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Abstract : We present a simple yet powerful technique for forming iterative methods 
of various convergence orders. Methods of various convergence orders (four, six, eight 
and ten) are formed through a modest modification of the classical Newton method. 
The technique can be easily implemented in existing software packages as suggested by 
the presented C"*"^ algorithm. Finally some problems are solved through the proposed 
algorithm. 



1 Introduction 

The most common and probably the most used method method finding a simple 
root 7, i.e. /{"j) = 0, of a nonlinear scalar equation 



fi^) = 0, 

is the Newton method. The classical Newton method is given as 



Xn+l 



n = 0, 1,2,3, . 
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It is well documented and well known that the Newton's method converges 
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251 : |26| . and references therein] . There exists many modifications 
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of the Newton method to improve convergence order 
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modifications of the Newton's method free of second or higher derivatives are 
activelyresearched. For example, third order convergent methods are presented 
in [ll; Il8; 20; 21], fourth order convergent methods are developed in [J; [1; [1; 
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[1 i; S 0; H [3 sixth order methods are developed in [H [li [13] and 
eight order methods are presented in [see|3l|, and references therein]. As there 
exists various modifications of the Newtons method. And from practical point 
(implementing these methods into a software package) one of the drawbacks of 
these wonderful methods is their independent nature. For example, if one has 
a software package which solves nonlinear equations by the well known fourth 
order convergent Jarrat method [1^. One may find it difficult to modify the 
existing software package to implement sixth order methods 28; [s^ and 
eight order methods [31]. 

In this work, we develop a technique. Which improves the order of conver- 
gence of the Newton method ^ from 2 to 2 x m. Here, m = 1, 2, 3, . . .. Thus 
through our scheme, one may develop 4th order, 6th order, 8th order, . . . con- 
vergent iterative methods. One of the beautiful fact of our scheme is that one 
needs modest modifications in the classical iterative method ([2]) for achieving 
higher convergence rates. And, which may be very effective when one wants 
to modify an existing software package for achieving higher convergence order. 
Let us now develop our scheme. 



2 The technique and convergence order of its 
various methods 

Before presenting our technique, first we will develop iterative methods of var- 
ious convergence orders. Let us consider the following 4th order convergent 
iterative method 



yn — 



liVn) 
f(Xn) 



1 + 2 



f{Xn) 



error equation for the above method is given as 
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Here, Ck = f^{l)/^-k and e„ = a;„ — 7. A proof of convergence of the above 
fourth order method is presented next. 

proof 

Using Taylor series of f{x) around the solution 7, and taking into account 
f{l) = 0, we get 



(6) 
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furthermore from the equation we have 

oo 



(7) 



fc=i 



and through a simple calculation we arrive at 

f{x) _ _ C2 2 _ 9 C3C1 - C2^ 3 _ 3C4C1^ - 7 C2C3C1 + 4 C2^ ^ 4 

J (a;) ci ci^ 



e„* + O (e„5) 
(8) 



Substituting ([5]) in (O yields 



C2 2 , r,C3Ci+C2 3 , 3 C4C1 + 7 C2C3C1 + 4 C2^ i , ( 5\ 

yn-l=—en +2 5 e„ H 5 e„ + O (e„ j. 
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Expanding /(?/«) around the solution 7 and using fO]), we obtain 
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From equations ([6|) and (fTO|) . we get 

/(;/«) _C2 2 C3C1 - 3 C2^ 2 _ -3c4Ci^ + IOC2C3C1 - 8C2' 
/(a;„) Cl ci^ ci-^ 



Now from equations ([8]), pT|) and ([4]), we find that 
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1 C2(-60C2^ + 12C 3C1)^^ ^ ^^^5, 

12 ci^ 



< + 0{el). 



(10) 



+ O (e„4) . 
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This proofs that the method ^ converges quartically. 

Let us now consider the following three step sixth order convergent iterative 
method 
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error equations for the above sixth order method is given as 



J_ C2 (-792 C3CiC2^ + 2160 €2'^ + 72 cg^ci^ 
72 
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Convergence order of the above method can be easily estabUshed through the 
Maple software package. We may notice that the method ([13]) requires evalua- 
tions of only three functions and one derivative during each iterative step. Let 
us now further consider the following eight order convergent iterative step 
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asymptotic error equation for the eight order method is given as 

1 C2 (77760 C2^ - 41472 cscicz" + 7344 c^^ci^Ci^ - 432 cg^ci^) 



(14) 
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We may notice that the eight order method ((T4)) requires evaluations of only 
four functions and one derivative during each iterative step. Based upon the 
similarity in methods ([4]), (fT3|) and ([Till . Let us consider the following method 

fi.Xn) 



f{x„) 



Pn 



Xn+1 
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through the Maple we verified that the above method is 10th order convergent, 
and error equation for it is given as 



C2 



2592 c? 



(2799360 C2* 



1959552 C3CiC2^ + 513216 cg^ci^cz'' 



-59616 cg^ci^ca^ 



2592 ei°. 



We may notice that the above tenth order method (|15|) requires evaluations of 
only five functions and one derivative per iterative step. 

Based upon the methods ([4]), p3)) . (fT4)l and (fT5|) . we conjuncture the exis- 
tence of the following scheme for generating iterative method of order 2 x m 
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2/2 
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It may be notice that a 2 x m order method formed by the above scheme will 
require m functions and one deravitive evaluation per iterative step. A C^^ 
implementation of the above scheme (I16p is presented in the Listing [TJ 

3 Numerical work 

The order of convergence ^ of an iterative method is defined as follows 



, ./(2/i) 

7(^n) 



lim 



(17) 



Here, e„ is the error after n iterations of the method. Then the approximate 
value of the computational order of convergence (COC) p [27[ can be find using 
the formula 

_ LogKx^+i -j)/iXn -7)1 

Log|(a;„ -7)/(a;„_i -7)!' 
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Listing 1: CH — h implementation. 

int main(){ 

mp ::mp_init(2005); 
mp_real tol = l.Oe— 300; 
unsigned int w = 5; 
unsigned int maxitr — 1000; 
std : : vector<mp_real> x(maxitr ,0.0) ; 
x[0] = -0.5; 
mp_real errl = 100.0; 
mp_real err2 = 100.0; 
mp_real gamma = —1.0; 
unsigned int m = 1; 
unsigned i = 0; 
unsigned int conv — 0; 
while(errl > tol || err2 > tol){ 
i f ( i >= maxitr ) { 
break ; 

} 

x[i+l] = x[i] — function (x [ i ])/ firstderivative (x [ i ]) ; 
mp_real xl = x[i+l]; 
mp_real xO = x [ i ] ; 

for (unsigned int k = 1 ; k < m ; -H-k){ 
mp_real y — x[i+l]; 

x[i+l] = x[i+l] - function (x [ i +1]) 
/ firstderivative (xO)*(1.0 

+ 2.0 * function (xl )/ function (xO ) ) 

errl = abs (y— x [ i + 1 ] ) ; 

err2 = abs ( function (x [ i + 1] )) ; 

if(errl < tol & err2 < tol){ 

std : : cout « "errl — " « scientific 
« setw(2*w) « errl 
« " err2 = " «setw(2*w) 
« err2 « std : : endl ; 

conv — 1 ; 
break ; 

} 

} 

i f ( conv ) break ; 

errl = abs(x[i+l] — x[i]); 

err2 = abs ( function (x [ i + 1] )) ; 

std :: cout « setw (w) « "itr. = " « i « scientific 
« " " « setw(2*w) « x[i+l] 
« std : : endl ; 

++i ; 

} 

mp_real rho ; 

for (unsigned int n = 1; n < i; -H-n){ 

rho — log(abs( (x[n + l] — gamma) /( x [ n] —gamma) ) ) 

/log(abs( (x[n] — gamma) /( x [ n — 1]— gamma ))) ; 
std :: cout « "rho = " « rho « std :: endl ; 

} 

mp : : mp_finalize ( ) ; 6 
return 1; 



All the computations reported here are done in the programming language 
C"*""*". For numerical precision, we are using ARPREC[32l. The ARPREC pack- 
age supports arbitrarily high level of numeric precision 32] . In the program (see 
the Listing 1) the precision in decimal digits is set at 2005 with the command 
mp::mp init(2005)|32l|. We have implemented the presented technique in the 
C"'""'' language. Listing [T] presents the main part of our implementation. 

For convergence, it is required that the distance of two consecutive approx- 
imations (|a;„+i — Xn\) be less than e. And, the absolute value of the function 
also referred to as residual be less than e. Apart from the conver- 
gence criteria, our algorithm also uses maximum allowed iterations as stopping 
criterion. Thus our algorithm stops if (i) \xn+i — Xn\ < e (ii) |/(x„)| < e (iii) 
itr > maxitr. Here, e — 1 x 10^^°°, itr is the iteration counter for the algorithm 
and maxitr — 100. See the C+"'" algorithm presented in the Listing [1] The 
algorithm is tested for the following functions [lH 



h{x) 


=x^ + Ax'^ - 


15, 




7 - 


i 1.347, 


f2{x) 


= sina; — x/3, 






7 - 


i 2.278, 


h{x) 


^lOxe-''^ - 1, 






7 - 


i 1.679, 


h{x) 


= cos X — X, 






7 - 


i 0.739, 


h{x) 








7 - 


i-1.000. 


k{x) 


^e^^ + cosx, 






7 - 


i 1.746, 


h{x) 


=Log(x2 +X + 2) 


-x-\ 


-1, 


7 - 


i 4.152, 


h{x) 


=arcsin(a;^ ^ 1) ^ 


x/2- 




7 - 


i 0.5948. 



We run the algorithm shown in the Listing[T]for four values of m: m = 1, 2, 3, 4. 
Here, in = 1 corresponds to the classical Newton method. We choose the same 
initial guess as found in the article (sij . So the reader may find it easier to com- 
pare performance of various methods. Table [1] reports outcome of our numerical 
work. The table [T] reports (iterations required, number of function evaluations 
needed, COC during second last iteration) for the Newton method (m = 1), 
fourth order iterative method (m = 2), sixth order iterative method (m = 3) 
and eight order iterative method (m = 4) . Computational order of convergence 
reported in the Table [1] was observed during the second last iteration. Following 
important observations were made during numerical experimentations 

1. In the Table [U the methods which require least number of function eval- 
uations for convergence are marked in bold. 

2. As reported in the Table[l] for the function f2{x), COC (during the second 
last iteration) is same for m = 1, 2, 3, 4. While COC for the second, third 
and fourth iterations is reported in the Table [2] 

3. From the Table [TJ we notice that for the functions f^ix), /4(x), frix) and 
fsix) the sixth order (w = 3) and eight (m = 4) methods requirs same 
number of iterative steps. Table [3] reports residual (|/(a;„)|) during the 
last iterative step. 
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/(■^) 


Xo 


NM(m = 1) 


TO = 2 


TO = 3 


TO = 4 


/l(^) 


1.6 


(9,18,2) 


(4,12,4) 


(2,8,5.66) 


(2,10,7.6) 


/2(X) 


2.0 


(23,46,1) 


(10,30,1) 


(7,21,1) 


(6,30,1) 


h{x) 


1.8 


(10,20,2) 


(4,12,3.99) 


(3,12,6.21) 


(3,15,8.22) 




1.0 


(9,18,2) 


(4,12,3.99) 


(3,12,5.90) 


(3,15,8.10) 


M^) 


-0.5 


(11,22,2) 


(5,15,3.99) 


(4,16,5.99) 


(3,15,6.75) 


k{x) 


2.0 


(9,18,2) 


(4,12,3.99) 


(3,12,5.99) 


(2,10.8.10) 




3.2 


(10,20,2) 


(4,12,3.99) 


(3,12,6.19) 


(3,15,8.19) 




1.0 


(10,20,2) 


(4,12,4.01) 


(3,12,6.35) 


(3,15,8.36) 



Table 1: (iterations, number of function evaluations, COC) for the Newton 
method (to = 1), fourth order method (to = 2), sixth order method (to = 3) 
and eight order iterative method (to = 4). 



fix) 


m 


= 1 


m, 


= 2 


TO, 


= 3 


TO, 


= 4 




10 


-i'Ah 


10 


-872 


10 


-671 


10" 


-lb'2-2 



Table 2: Computational order of convergence (COC) at the second iteration. 



fix) 


m = 1 TO = 2 TO = 3 TO = 4 


f7{x) 

fs{x) 


IQ-^T^ 10"^'''^ IQ-IS'2'2 
10-347 1 0-744 10-800 10-1302 



Table 3: Residual (|/(a;„)|) 
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